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A three-body calculation of the reaction 'yd rjnp in the energy region from threshold up to 30 
MeV above has been performed. The primary goal of this study is to assess the importance of the 
three-body aspects in the hadronic sector of this reaction. Results are presented for the 77-meson 
. spectrum as well as for the total cross section. The three-body results differ significantly from those 

' predicted by a simple rescattering model in which only first-order rjN- and A^'A^'-interactions in the 

final state are considered. The major features of the experimental data are well reproduced although 
right at threshold the rather large total cross section could not be explained. 
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I. INTRODUCTION 

In this letter we would like to present theoretical results for total and differential cross sections of the reaction 
jd — > r]np close to the threshold region. The need for an improved theory for describing the low-energy ry-production 
I on the lightest nuclei stems from the recent precise measurements of such processes on light nuclei using hadronic 
and electromagnetic probes. The experimental results are substantially larger than a mere phase-space calculation 
, near threshold It is reasonable to associate this effect with a rather strong attractive force which governs the 

' dynamical features of the ryA^ A^-system in the low-energy regime. In a previous paper we found that the rjNN- 
(-H , dynamics allows the existence of virtual s-wave three-body states near zero energy. The natural extension of these 
findings would be to investigate to what degree such states may influence the processes considering the TyA^A^-system 
in different spin-isospin channels. Of particular interest are the reactions of coherent (7c? rjd) and incoherent 
(7^ — > r]np) 77-photoproduction on the deuteron, which recently have become the subject of attention of different 
experimental collaborations |^,^,^. 

As for the coherent process, the investigation of this channel within a three-body approach has been performed in 
1^ . Unfortunately, due to the smallness of the coherent 77-photoproduction cross section resulting in extremely small 
77-meson yields, only qualitative statements about the size of the experimental cross section of the coherent reaction 
5J] ' is available at present. 

I In this work, we focus our attention exclusively on the incoherent 77-photoproduction where the deuteron breaks up. 
Recent precision measurements of the break-up channel require a more sophisticated theoretical analysis. Indeed, 
this reaction has already been considered in a previous paper within the so-called rescattering model as a first 
step beyond the simple impulse approximation (lA). In this model, final state interaction is taken into account by 
including the single rescatterings between two of the final particles as most important correction to the lA. As was 
shown in Q], without such rescattering effects a strong suppression of jd — > rjnp cross section was found because 
the lA requires a large spectator nucleon momenta for producing a low-energy 77-meson on the deuteron. For this 
reason, the inclusion of first order rjN- and AAf-interactions provides a mechanism to distribute the high momentum 
between two nucleons and thus avoids this suppression. In view of the importance of the first-order rescattering it is 
natural to ask about the role of higher order rescattering terms. Indeed, in a system of three low-energy particles with 
rather strong attractive forces, one can expect a considerable overlap of their wave functions. Therefore, there is no 
guarantee that the contributions from higher order rjN- and A^A^-rescatterings will be small. This fact distinguishes 
the present reaction from those with pions where the smallness of the pion-nucleon scattering length allows one to take 
into account only the lowest order A^A^-interaction in the final state just above threshold g. Thus, the rescattering 
approximation can hardly be considered as a satisfactory method for describing low-energy 77-photoproduction on light 
nuclei. It must also be noted that the successive inclusion of the next higher order terms of the multiple scattering 
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expansion is not well advised, because the corresponding Neumann series is expected to converge very slowly due to 
the proximity of the virtual three-body poles mentioned above. As a consequence, the inclusion of terms of quite 
high order will be required rendering such an approach inadequate. These considerations point to the necessity of a 
three-body approach, where the final state interaction is included to all orders. 

In the present paper, the reaction 7c? — > rjnp is considered as a three-body problem in the hadronic sector. The 
formalism is based on a separable s-wave approximation for the two-body potentials. This approach may be justified 
by the fact that in the energy region under consideration the yyiV-resonance 5*11(1535) as well as the ^6*0 and ^Si poles 
of the NN scattering matrix determine essentially the driving r]N- and iViV-forces. Our primary goal is to investigate 
the main properties of the incoherent ry-photoproduction on the deuteron near threshold. We also compare our results 
with those predicted by the first order rescattering approximation in order to assess the shortcomings of this method. 

In the Sect. |l| the three-body equations for the rjnp transition amplitude are derived. In Sect. [II 



calculate the meson spectrum for a fixed photon energy and compare our results for the total cross section with 
available experimental data. 



II. GENERAL FORMALISM 



Here we will briefly describe the main points of the formal part, and comment on the most important approximations 
involved in our calculation. A more detailed and rigorous development of the formalism as well as the method for 
the numerical solution of the three-particle integral equations will be presented elsewhere. Firstly, we note that the 
three-body concept of the 7y-photoproduction on a deuteron possesses in principle a straightforward solution of the 
three-body multichannel problem where the coupling between the channels jNN and r/NN is taken into account. 
However, since we will restrict ourselves to the first order only in the electromagnetic interaction, it is more appropriate 
to obtain the dynamical equations of our model starting from the Faddeev formalism for the pure hadronic rjNN 
system. Our approach with respect to the hadronic part is essentially based on the separable approximation scheme 
of Alt, Grassberger and Sandhas |^ with an isobar ansatz for the rjN- and iViV-interactions. As is mentioned above, 
this ansatz should represent a fair approximation of the real dynamics in the energy region examined here. We use 
the following notations for the basic three-body channels: 

N* for the 7?A''-isobar plus spectator nucleon, 

d for the TV TV-isobar plus spectator meson. ^ ' 

and analogous notations for the two-body subsystems. Then the separable input 

U^\h)r,{M {t^N*,d) (2) 
allows one to obtain two coupled operator equations for the rearrangement amplitudes 

Xd = 2ZdN*T]y.XN* , (3) 
Xn* — Zjyd + ZN-dTdXd + Zn'N'TN'Xn' ■ 



The amplitudes X^ and X^f determine the transitions from the initial on-energy-shell state rjd to the quasi- two-body 
final states r/c? and NN*, respectively. The potentials Zij have a conventional meaning, namely 

Z^, = iMGolf,) , (4) 

with Go being the free three-body Greens function. 

Turning now to the photoproduction process jd — > rjnp, we treat the electromagnetic interaction perturbatively, 
keeping only the terms up to first order in the e.m. coupling. Therefore, in order to introduce the electromagnetic 
part, one just has to replace in (||) the ?7-meson by the photon in the entrance channel. This yields then the coupled 
equations which are formally similar to 

Xd = 2ZdN*TN*XM' , (5a) 

Xn" = B]y. + Zjy.dTdXd + Zn'N'TN'Xn. . (5b) 



Here the "born term" Bn" is determined as 



Br,^={/^}\G^^^\fd), (6) 
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where G^^^ is the free propagator in the 77VA^-sector and 1/^2) is the electromagnetic vertex function for the transition 
— !■ A^*. Equations are represented diagrammatically in Fig. |l|. Clearly, the impulse approximation is equivalent 
to the replacements Xd and Xn* — > in the right-hand sides of (||), i.e., 

Xl^ = 0, (7a) 
Xj^, = . (7b) 

In order to obtain the first-order rescattering approximation, in which only the single rjN- and A^A^-scattering terms 
are retained, one has just to set Xd — and Xn* = Bn* ■ Then one has 

j^rcsc _ 2ZdN'TN'BN* , (8a) 

X]^l'^ = Bn' + Zn,^,tn'Bn, . (8b) 

After projecting (||) and correspondingly (|^) and (^ onto the L = states, we end up with a system of one- 
dimensional integral equations in momentum space, which can be solved by matrix inversion. The well known 
problem of logarithmic singularities by the rearrangement terms ZjvAf and Zj^td was avoided using the contour 
deformation method of . 

The transition matrix element for the break-up process jd rjup may be expressed in terms of the rearrangement 
amplitudes (||) as follows 

{pi,P2,q\TiW)\k} = {p„i\fN'} TN'iW-Ep, - {p2\XN'{W)\k) 

2 

+ {p,2\fN') TN'iW-E,, - ^^) {pi\XN^(W)\k) (9) 
+ (Pl2|/d) TdiW -LO^- ^) {q\Xd{W)\k) . 

Here, W denotes the total three-body cm. energy and k the photon momentum. Furthermore, the three-momenta of 
the final nucleons and ry- meson are denoted by pi, p2 and q, and their total energies by Ep^ and w^, respectively. 
The propagators TN*{E]\[t) and Td{Ed) depend on the invariant energies of the corresponding two-body subsystems. 
Their explicit form is given in The arguments of the vertex functions {p\fi) are the respective two-body relative 
momenta, whose nonrelativistic expressions read {i — 1,2) 

qMN-Ptfrin ^ Pi - P2 
" ~U — \ ' ^'12 = 7, • 10) 

In the actual calculations, only the s-wave configuration ( J'^; T) = (0~; 1) of the yyA'^iV-system was taken into account. 
As was discussed in [Q, due to the spin-isospin selection rules, this state gives the major contribution to the —> rjnp 
cross section. We also use the nonrelativistic kinetic energies for all three final particles, since we restrict ourselves 
to the near-threshold region where this approximation is well justified. The contribution from pion exchange in the 
rearrangement potential Zjvat* was ignored in view of its insignificance as was shown in ||^. The parametrization of 
the hadronic separable t-matrices (||) is given in detail in |^ . For the s-wave singlet and triplet A^-states we use the 
fit of Yamaguchi |jll|. The A^* parameters were chosen in such a way that the main hadronic width of the 5ii(1535) 

resonance, Tn^^t^n ~ ^N^-triN ~ 75MeV [^, is reproduced. As for the electromagnetic vertex l/^i), it is natural 
to associate it with the 5*11(1535) helicity amplitude A^fJ^ which is of great experimental importance. Thus 

one has 



Since we restrict our consideration to the rjN N-sta.te with isospin T = 1, only the dominant isovector part of the 
7A'' 511(1535) transition has been taken into account. The corresponding amplitude was derived as follows. First 
we have determined the helicity amplitude in the proton channel ^^72 fitting the experimental cross section data 



for the elementary process ^p —>■ rjp of \13\. Then the isovector part ^[^2 extracted as ^^^2 ~ 0-9^i/2 which is 



compatible with different analyses [Qi^jj^]- 

Finally, we present here the formula for the exclusive cm. cross section used in the present calculation 
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^ (2^)5 '""V 2 ^"'^'^^1' ^'^"^ ^^^^ ' ^^^^ 

where i?d is the deuteron energy and is the azimuthal angle of the nucleon momentum pi counted from the k-q 
plane in the frame with the z-axis directed along q. The factor i appears after summing and averaging over the spin 
projections. The semi- inclusive cross section da'^ ^ /iiw^ and the total cross section atot considered below are obtained 
from (|l^) by appropriate integrations. 



III. RESULTS AND DISCUSSION 



We would like to beginn our discussion with the ry-meson spectrum da'^'™' / duo^ as function of the rj cm. energy, 
shown in Fig. ^ for two different photon energies. Its specific form was discussed already in ||^. A strong s-wave 
attraction between the two nucleons at very low relative energy shifts visibly a major part of the spectrum towards 
higher energies, resulting in a rather pronounced peak near the boundary of the phase space at the maximum of the 
available meson energy, where the relative energy of the two nucleons is close to zero. It is obviously a manifestation 
of the large NN singlet scattering length, i.e., of the antibound So-state. As our calculation shows, the first-order 
rescattering approximation underestimates substantially the complete three-body result just above threshold, but 
overestimates it at higher energies. Although the form of the spectrum is not changed qualitatively, the peak is much 
more pronounced in the complete calculation. 

The difference between the results of the first-order rescattering and the complete three-body calculation is most 
apparent in the energy dependence of the total cross section shown in Fig. |^, where we compare it also with the 
inclusive ^d rjX data from |^ . As shown in ||^ and noted in the introduction, rescattering gives a very important 
contribution in particular just above threshold, where it leads to a strong enhancement of the meson yield. However, 
truncating the multiple scattering expansion after the first order gives only a qualitative description of this effect. 
Inclusion of rescattering to all orders within the three-body approach results in a considerably stronger effect in the 
near threshold region so that the energy dependence of the cross section becomes slightly concave (solid curve "1" 
in Fig. which is not characteristic for photoproduction reactions with more than two particles in the final state. 
We associate this effect with the virtual r]NN-pole found in |^ in the {J^;T) = (0~; l)-state which pulls down the 
essential part of the meson yield to the near-threshold region. 

When comparing our results with the experimental data we see, first of all, that the three-body calculation consti- 
tutes a very important improvement of the theory which reproduces fairly well the qualitative features of the data in 
the near threshold region up to about 650 MeV photon energies, in particular with respect to the energy dependence. 
On the other hand, one notes right at threshold still a substantial underestimation of the data by the theory. One 
reason for this disagreement might be, of course, a contribution of the coherent channel 7c? rjd to the experimental 
inclusive yield. However, to ascribes the whole difference to the coherent process alone would mean to assign a rather 
large cross section of about 100 nb at £'^=635 MeV to the coherent process. While such a value is not in conflict 
with the recent measurements of [|[ it nevertheless overestimates substantially the corresponding PWIA-calculations 
for the case = El- In order to reach a more definite conclusion one needs, of course, more precise 

data for the total cross section of the coherent reaction. Finally, with respect to the underestimation of the data at 
energies above 650 MeV, one should keep in mind that we have included in the calculation only the j]NN-state with 
total orbital momentum L^O. Thus at higher photon energies, where the higher partial waves are needed to fill the 
available phase space, our calculation cannot be considered to be realistic. 

In order to complete our discussion, we have investigated the sensitivity of the results to the strength of the 
interaction, i.e., to a variation of the 7V*-parameters. As is shown in Q the rjNN* and ttNN* coupling constants 
influence the position of the r/NN S'-matrix poles rather strongly and thus may also affect the cross section value. 
Indeed, this is proven by a calculation using different sets of the rjNN*- and TrA^iV* -couplings, the results of which 
are also shown in Fig ^ In order to illustrate the size of the variation of the interaction strength by varying the model 
parameters, we quote also the corresponding values of the TyA^-scattering length arjN associated with the choice of the 
coupling constants. This parameter characterizes the low-energy ryA'^-interaction and is quoted in most analyses of the 
TyA^-interaction (for a review of the a^Ar-values see, e.g. [|5j). We conclude from Fig || that the results are not very 
sensitive to the choice of a,,Ar, so that even with a rather large value of Qt^n = (0.91 -I- i0.25)fm (solid curve "2" in 
Fig. H), we cannot describe the data sufficiently well. 
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IV. CONCLUSION 



We have presented theoretical results for the cross section of the reaction rjup in the energy region up to 30 MeV 

above threshold. The calculation is performed within a three-body formalism with s-wave separable rjN- and NN- 
potentials. We found that the first-order rescattering approximation is insufficient, since the full three-body treatment, 
including the full hadronic final state interaction, enhances significantly at lower energies the ry-meson spectrum while 
at the higher energies a less pronounced but still sizeable decrease appears. These results demonstrate clearly that 
any " minimal" realistic model for the 77-production near threshold must be based on a three-body description of the 
77 A^-dynamics . 

Our calculation explains at least qualitatively the anomalous behaviour of the experimental 7y-meson yield observed 
just above threshold. On the other hand, a quantitative agreement with the experimental results is not achieved, 
even with a rather large r^iV-scattering length a,jAr=(0. 91-1-10.25) lying on the boundary of values allowed by modern 
analyses. Whether the discrepancy could be explained by the contribution from the coherent channel 7c? — *■ rjd is at 
present an open question. 

In view of our results on the 7/-meson spectrum as presented in Fig. ||, we would like to emphasize the fact that 
very useful and interesting information on the mechanism of the low-energy 77-photoproduction may be obtained 
from a measurement of this spectrum. The characteristic feature of the meson spectrum with its sharp, pronounced 
maximum close to the high energy kinematic limit should make it possible to identify the final state interaction as a 
major source of an anomalously large cross section of the reaction rjnp close to the threshold. On the contrary, 

if there is no simple explanation for the noted shift between the data and the model predictions, then one has to 
invoke some additional unexplored mechanism in 77-photoproduction. 
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FIG. 1. Diagrammatic representation of the three-body equations for ry-photoproduction on the deuteron. 
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FIG. 2. The r;-meson spectrum for the reaction -yd — > -qup versus the total cm. energy cu^ for two lab photon energies. 
Notation of the curves: dotted: impulse approximation (lA), dashed: inclusion of first-order rjN- and A''A'^-rescattering, solid: 
full three-body calculation. 
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FIG. 3. Total cross section for the reaction yd rjup. Notation of the curves as in Fig. Q. The results obtained within 
the three-body approach with different sets of rjNN* and ttNN* couplings are presented as the solid curves "1" and "2". The 
corresponding values of the 7;A''-scattering length ar,N are given in the legend. The inclusive ■yd — > r/X data are taken from Q. 
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